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In order to more effectively contribute to the 
research program of the Instrument Physics Research Section, 
NASA Langley Research Center, work under this grant was 
directed toward the problem of measuring diffusion para- 
meters for the beryllium impurity in silicon. The very 
-5 2-1 

large (~ 10 cm sec ) diffusion coefficients encountered 
in this study required the determination of impurity pro- 
files for thick samples (~ 1 cm), due to the deep penetration 
of the impurity occurring during reasonable diffusion times. 
Since conductivity values are readily converted into con- 
centrations of electrically active impurities, the major 
problem became that of accurately determining the conduc- 
tivity profiles of beryllium diffused silicon samples. 

Details of the problem and its solution are contained in 
the attached report. 
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Profiles of Depth Dependent Conductivity 
As Determined by Electrical Surface Probes' , 

R. K. Franks 

Physics Department, Appalachian State University, Boone, North Carolina 

AND 

J. B. Robertson 

NASA,- Langley Research Center, Hampton, Virginia 

Four-point probe measurements on samples having depth 
dependent conductivities are interpreted in terms of the 
conductivity profiles. An exact solution of the problem of 
exponentially depth dependent conductivity serves as a basis 
for this interpretation. Applications include surface conductivity 
determination where the form of the conductivity profile is known, 
and conductivity profile determination from probe measurements 
taken as the sample surface is progressively lapped away. 
Application is limited to samples having conductivity monotonically 
decreasing with depth fran the probed surface. 
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INTRODUCTION 

Electrical surface probes have proven useful in geological 
1 

exploration as well as in conductivity profile determinations of 

2 3 

semiconducting samples . ’ Published methods of interpreting surface 

probe measurements in terms of a depth dependent conductivity are based 

upon a layered model with conductivity values constant within each layer. 

These methods use solutions of Laplace's equation within each layer with 

the implicit assumption that charge accumulations at the boundary are 

excluded from the layer volume. This neglect of the physical requirement 

that the charge accumulation occupy a non-zero sample volume gives rise 

3 

to an error of undetermined magnitude. A second disadvantage of these 

methods is their complexity which makes consideration of more than a 

few layers prohibitively difficult. 

The four-point probe as illustrated in Fig. 1 is often used in 

measurements of conductivity. Among its many advantages are ease. of 

application and, for homogeneous or very thin samples, ready inter-' 

U d 

pretation. ’ The problem of a slab of finite thickness with an 
exponentially depth dependent conductivity has been recently solved,^ 
and this exact solution serves here as a basis for interpreting four- 
point. probe measurements on samples having other forms of depth dependent 
conductivity. Applications include surface conductivity determination 
where the form of the conductivity profile is known , and conductivity 
profile determination from four-point probe measurements taken as the 
sample surface is progressively lapped away. In this paper we present 
a simple method of conductivity profile analysis which appears to be 
reasonably accurate for samples having a conductivity monotonically 
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decreasing with depth from the probed sample surface. 

In the following treatment, some of the more important assumptions 
and stipulations are: 

(1) The probe is made up of four coll inear contacts with point 
spacing S. 

(2) The two inner contacts are used to determine the potential 
difference, V (volts), due to the current, I (amperes) > through the 
outer pair of contacts . 

(3) Sample conductivity, <r(z) £ohm-cnr| , is a function only 
of depth below the probed surface. 

(U) The probe point spacing is small compared to the probe distance 
from the nearest lateral boundary. 

EXPONENTIALLY DEPTH DEPENDENT CONDUCTIVITY 

6 

We have shown previously that for conductivities, a, of the form 

tf(z) = <x(0) exp(-az) , 0) 

the probe potential difference may be expressed as 

V =[l/2TrSCj(0)] F aS (T) , (2) 

where z is the depth below the probed surface, T is the sample thick- 
ness, and a is a constant characterizing a particular exponential depth 
dependence of the conductivity. F (T) is an appropriate correction 

3.0 
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factor given as 


v O'* - aS jgU 

J aS U; " 1-exp (-aT) 


Us 


Jjn^Tr^ 


n4 VUn 2 tr 2 +a-T 2j 



where the are Bessel functions of the zeroth order. 

The usefulness of the preceeding relationships is greatly enhanced 

by a recently pseudo-derived approximation to F (T) which is good to 

3.0 

within 0.1$ error for all a 20. Our approximation is 


•p, ( T w_aSjnU_ 

aS K J ~ 1-exp (-aT) 

+ ^r-CST'^nU) - la] S J?n z] Tanh/ ( /a| S) 1 [exp( |a/ T)-1 2> 
|aj S jexp( |aj T) -ijf 1 ^nU + exp(-ST iUiU) 

which for the case a=0 simplifies to 


F^tTj^ST- 1 ^nU ♦ e^ST- 1 SnU) Taah (S^T) 
* ST 'JZnU + exp(-ST j?nU) 


(U) 


( 5 ) 


The restriction of the approximation to a'Z'O may be appreciated by 
verifying that the leading terms of Eqs. (U) and (3) are identical so 
that the second term of Eq. (U) serves as an approximation to the 
summation term of Eq. (3)- The approximated second term deteriorates 
as Ja{ increases. However, for positive a, this deterioration is 
compensated by the increasing dominance of the first term. Unfortunately, 
for negative a the approximated term becomes increasingly dominant as 
N increases. 
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The function F _(T) decreasestoward. a limiting, value as the 
ab 

thickness, T, increases. This implies that the degree of involvement 
of the sample material in the conduction process decreases with in- 
creasing depth from the probed surface. For sample thicknesses small 
compared to the probe spacing, the effect of sample material upon the 
probe potential is 'almost independent of depth below the probed surface 

and the leading term of F (T) dominates. In this case ,Eq. (2) reduces 

ab 

to the expression^ 

V = (j?n2/ rr) (3/cr ) , (6) 

"S 

appropriate for a thin infinite sheet having sheet conductivity, ^ . 

For our purposes, we find it convenient to extend the range of Eq. (6) 
to any thickness by introducing a multiplicative factor, W (z), into 

clb 

the usual definition of sheet conductivity so that 

2 s ' o/ Tw a S (z) S(z)dz - (7) 

Combining Eqs. (6) and (7) to obtain 

W aS (z)cr(z)dz = (in2/n) (I/V) , (8) 

we observe that a knoitfledge of (I/V) as a function of either S or T 
uniquely determines vr(z) . 

Since Eqs. (2) and (6) both describe the same situation, we must 


have 
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o/ T W as (z)rr( z )dz = S M j(0) F a g(T) . (9) 

We assume that the degree of involvement of the material at a given depth 
in the conduction process is solely dependent upon the conductivity and 
thickness of the overlying material. The mathematical implication of 
this assumption is that the weighting factor W (z) , in the left hand 

3.0 

term of Eq. (9) is independent of the limit of integration, T. Thus 
we may differentiate the terms of Eq. (9) with respect to T to obtain 

W aS (z) = S £nkg[(0)a~\z) jdF~g(z)/dzj . (10) 

Some representative weighting functions appropriate to samples having 

exponentially depth dependent conductivities are shown in Fig. 2. 

While the assumption of the preceeding paragraph is not strictly 

valid, it appears to be a good approximation to reality in the case of 

semiconductor slabs having the conductivity decreasing as the depth 

below the probed surface increases. It may be helpful to examine an 

equivalent assumption that the conduction process at depths greater 

than z has no effect upon the electric field at depth z. This clearly 

suggest that our assumption is far from valid when the underlying 

material is of relatively high conductivity. 

For interpreting four-point probe measurements made as the top 

sample surface is progressively lapped away, we are pleased to note 

that W„ (z) serves as a useful approximation to all W (z), a^O. 

Uo 

Replacing w a g( z ) by W^( z) in application of Eq. (8) to constructs 
devised using a wide range of exponential dependencies, the maximum 
error encountered in surface conductivity, cr(0), was less than 6%. 
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This magnitude of cumulative error has little effect upon the form of 
the profile . The maximum error occurs for large sample thickness . and 

-1 

at a conductivity profile characterized by . the exponential factor a=S 

This error in sample surface conductivity goes rapidly to zero as a 

-1 

approaches zero and has decreased to 0.6 % at a=5S . This is consistent 
with the observation that the true weighting function, ¥ (z), will be 

clO 

significantly higher than ¥ (z) only when the material overlying the 

vO ' 

region about z has a relatively high conductivity; or equivalently 
only when the material about z has a relatively low conductivity so that 
the effect of the material at z upon the probe potential is relatively 
unimportant . 

APPLICATION 

The modified concept of sheet conductivity, as seen in Eq. (7), 
can be readily extended to other forms of depth dependence of conductivity, 
subject only to the availability of an appropriate weighting factor. 

For samples not having exponentially depth dependent conductivities, the 
weighting function, ¥ (z), must of necessity be an approximation. The 

O 

use of the weighting function ¥ (z) for samples having conductivity 

profiles characterized by a introduces no computational error and one 
might expect that for samples having conductivity profiles similar in 
form to an exponential depth dependence, the use of an appropriate 
¥ (z) would produce highly accurate results . 

cLq 

The primary motivation for this study was the need of a convenient 
and reliable method of determining concentration profiles of impurities 
diffused into semiconductors. Conductivity values are readily converted 
into concentration's of electrically active impurities using mobility 
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data such as that summarized by E. M. Conwell. In silicon, typical 

diffusion depths of substitutional impurities are a few microns and the 

easily applied four-point probe method offers high resolution and ready 

£ 8 

interpretation. ? Pentration depths of the fast diffusing interstitial 

impurities are measured in millimeters, however, and require use of a 

more complex interpretation of four-point probe data. 

To accomplish a conductivity profile analysis we may consider the ,, / 

/ ] 

sample to be an N-layered structure as illustrated in Fig. 3. Layers 
are characterized by an average conductivity _cFj_, a thickness (AT) i? and 
an average weighting factor W The generalized form of Eq. (8), in 
a form suitable for numerical integration, is 

N _ 

J w si % (AT) i = (j?n2/TT) (I/V) s * (11) 

By applying Eq. (11) to surface measurements taken with N different 
probe-point spacings, we obtain N equations with N unknown 

Calculation is simplified in cases where four-point probe measure- 
ments are made on the top sample surface at each . stage as the layers 

are progressively lapped away. In this case we use a constant probe 

"btl 

spacing and obtain the average conductivity of the j layer as 

3 j * jW/n) (X/V) ^W^i. ( iT )l ^(41) J" 1 . . (12) 

It should be clear that the calculations must begin with the sample 
remnant and work back toward the first layer to be removed. 

It is sometimes desirable to find the value- of the surface conductivity 
of a sample for which the form of the depth dependent conductivity profile 
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is already known. Rewriting <j( z ) as g(0)f(z) in Eq. (8) and solving 
for the surface conductivity, cr(O) , we obtain 

<?(0) = (I/V) s (ln2/v)j/ 0 f T W g (z)f(z)dz . (13) 

t 

For application of Eqs. (12) and (13) we take ¥ (z), the 

U j 

weighting function appropriate to a homogeneous sample, as a fair • 
approximation to the true ¥ (z) for all samples having conductivities 

L> 

decreasing monotonically with increasing depth from the top surface. 

As in our previous discussion of the ¥ _,(z), this assumption is consistent 

with our expectaticn that the true ¥ (z) will be significantly higher 

o 

than ¥_ (z) only when the overlying material has a rather high conductivity 

v J 

relative to the material at depth z, so that the effect of the material 

at depth z upon the probe potential is relatively unimportant. 

In addition to the theoretical contructs involving constant and 

exponentially depth dependent conductivities, conductivity profiles due 

to gaussian and complimentary error function distributions of impurities 

in silicon have been employed to verify the usefulness of Eq. (12) in 

making profile analyses. Eq. (12) is being used extensively in studies 

of the beryllium impurity in silicon and yields impurity concentration 

profiles consistent with those expected from imposed experimental 
9 

conditions and with those obtained by mass spectroscopy analysis. 

In Fig. k we see the results of a 27 layer determination of the 
concentration profile of electrically active beryllium as a diffused 
impurity in silicon. The beryllium was introduced into the silicon at 
high temperature from one side only. The experimentally determined 
concentrations are shown fitted to the expected complimentary error 
function profile. 



10 


List of References and Footnotes 

Work supported by the National Aeronautics and Space Administration. 

■'p. K. Bhattacharya and H. P. Patra, Direct Current Geolectric 

Sounding (Elsevier Publishing Co., Amsterdam, 1968). 

^T. H. Yeh and K. H. Khokhani, J. Electrochem. Soc. , 116 . 1 U6l (1969). 

3 

P. A. Schumann, Jr. and E. E. Gardner, J. Electrochem. Soc., 116 , 

87 (1969). 

^L. B. Valdes, Proc. IRE, h2, U20 (1 95U) - 

'’f. M. Smits, Bell System Tech. J. , 37, 711 (1958). 

^R. K. Franks and J. B. Robertson, NASA Tech. Note. TN D-650k , Oct., 

1971. 

7 E. M. Cornell, Proc. IRE, U6, 1281 (1958). , ■ 

8 

E. Tannenbaum, Solid St. Electron., 2, 123 (1?6l). 

9 

Research Triangle Institute, Integrated Silicon Device Technology. 

Vol. IV — -Diffusion , ASD-TDR 63-316 Vol. IV, Durham, North. Carolina, 1961;. 



11 


Captions for Figures 

Fig. 1 . Four-point probe on flat semiconductor. 

Fig. 2. Some representative weighting functions, ¥ (z) , appropriate 

for samples having depth dependent conductivities of the form 
t(z) = o-(O) exp(-az). Units are chosen such that S = 1. 

Fig. 3. Multi-layered sample geometry assumed for conductivity profile 
determinations . 

Fig. U. A 27 layer determination of the. concentration profile, N(z), 
of electrically active beryllium as a diffused impurity in silicon. 












